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We show that the asymmetric tunnehng spectrum observed in the Cuprate superconductors stems 
from the existence of a competing order. The competition between the competing order and su- 
perconductivity can create a charge depletion region near the surface. The asymmetric response of 
the depletion region as the function of the external voltage causes the asymmetric tunneling spec- 
trum. The effect is very general in a system which is near the phase boundary of two competing 
states favoring different carrier densities. The asymmetry which has recently been observed in the 
point-contact spectroscopy of the heavy fermion superconductor CeCoIns is another example of this 
effect. 
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I. INTRODUCTION 

One important characteristic of the strongly correlated 
systems is the possibility of the existence of many com- 
peting orders. To identify these competing orders and 
study the interplay between them are two major tasks 
in this field. In high temperature cuprate superconduc- 
tors, competing orders, such as magneticfi*^ stripeji^'^ 
d-density wave^^ microscopic current orders and so 
on, have been proposed. Although it is still a matter of 
debate, the competing orders have been widely consid- 
ered to be responsible for the unusual properties and the 
universal phase diagrams of the materials. The competi- 
tion between superconductivity and antiferromagnetism 
or other orders is also present in other strongly correlated 
electron systems such as heavy fermions. 

In the cuprates, although the static global competing 
phases are rarely detected and the direct transition be- 
tween the superconducting and the competing orders has 
not been observed, recent STM experiments have sug- 
gested that a static order may exist iSiiSiii The compet- 
ing phases are expected to exist in the underdoped re- 
gion. In the multiple layer cuprate materials, the middle 
layer is naturally underdoped. In this case, the com- 
petition between the competing orders and superconduc- 
tivity has been used to explain the behavior of Tc as 
a function of the number of layersi^^ In this paper, we 
show that the carrier distribution close to the surface 
can also be naturally inhomogeneous along c-axis and 
the competing order is generically enhanced at the sur- 
face. The competition between the competing order and 
superconductivity can lead to an asymmetric tunneling 
spectrum. The effect could be the natural explanation of 
the asymmetry universally observed in the STM experi- 
ments in the cupratesl^ and the metal contact tunneling 
experiments^^ in heavy fermion systems. 

A particle-hole asymmetry is observed universally in 
different cuprate compounds. In traditional BCS type 
superconductors, tunneling spectrum is expected to be 
symmetric between negative and positive bias voltage. 
Such a large and universal asymmetry observed in high 



temperature superconductors remains a challenging puz- 
zle to be solved. A comprehensive theory to understand 
the phenomena may be tied to the ultimate understand- 
ing of high temperature superconductivity since there are 
only a few universal features identified crossing different 
types of cuprate materials. A few theoretical ideas based 
on the t-J model have been proposed liSiiLi^ 

In this paper, we propose that the origin of this 
asymmetry is indeed tied to one fundamental physics 
in cuprates, the existence of the hidden competing or- 
der. We show that the asymmetry stems from the com- 
peting order at surface, which has been detected by the 
same STM measurement >2ii2iii When an external voltage 
is applied, the competing order near the surface can be 
weakened or strengthened depending on the sign of volt- 
age. Therefore, the region of competing phases can be 
asymmetrically responded to the external field. This ef- 
fect leads to asymmetric tunneling matrix elements. Such 
an effect is expected to be universal as soon as the com- 
peting orders favor different doping concentrations from 
the superconducting states. Interestingly, we notice that 
the similar effect indeed has been observed in tunnel- 
ing experiments in heavy fermion systems^ The sample 
used in the experiment is also very close to the phase 
transition boundary between the AFM and SC orders. 

In a semiconductor, charges can be accumulated on 
the surface due to the electronic surface states. In the 
bulk close to the surface, the space charges can be in- 
duced to screen the surface charge. This effect creates a 
space charge layer at the semiconductor interface. In the 
cuprates, even without electronic surface states, a natural 
charge depletion region near the surface can be created 
since the carriers are introduced stoichiometrically. The 
number of the carriers induced in the Cu-0 layer close 
to the surface from the neighboring layer A-0 where A 
is La (Sr) in LSCO (BSCCO) can be very different from 
that inside the bulk. On the surface, free charges can 
be attached to the A-0 plane to reduce the electrostatic 
energy and effectively reduces the doping concentration. 
Such an effect leads to the lower carrier density in the 
Cu-0 plane close to the surface. In Figure ^ we plot 
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FIG. 1: Schematic picture of the carrier distribution and 
electrostatic field distribution along the c-axis in a double 
layer structure. The solid and darker solid lines represent 
the A-0 and Cu-O planes respectively and the dash arrows 
represent the electrostatic field, (a) and (b) represent the 
distributions in the bulk and close to the surface respectively. 



a schematic picture to show that in the presence of the 
surface, the electrostatic field can be reduced close to 
the surface to save the electrostatic energy vifhich is the 
origin of the doping inhomogeneity in the layered com- 
pounds of the type Tl2Ba2Ca„_iCu„02ri+4 and related 
structures. The effective carrier contribution from the 
A-0 layer on the surface is diminished in order to save 
the electrostatic energy. As a result, the doping concen- 
tration on the Cu-O layer close to the surface is lower 
than that in the bulk. Therefore, in the hole (electron) 
doped cuprates, a hole (electron) depletion region can 
be created near the surface. In the semiconductor, it is 
well known that an external voltage can change the space 
charge distribution near the surface and the depletion re- 
gion has asymmetric response to an external electric field. 
Thus one can expect that similar effect happens in the 
cuprates, even though it is not clear how the asymmetric 
response could create a large asymmetric STM tunneling 
spectrum if the sample is simply in the superconduct- 
ing state. In the presence of competing orders, however, 
the effect could be crucial. The variation of the carrier 
distribution could lead to a change of the configuration 
of the competing orders near the surface. In particular, 
when the competing orders are insulating, the tunneling 
amplitude can be significantly modified by the change 
of the configuration of the competing orders. Such an 
effect could be accounted for the observed asymmetric 
STM tunneling spectrum. 



II. MODEL AND FORMULATION 

We start to discuss about the physics in the struc- 
ture shown in Figure ^ The total charge carrier energy 
near the surface is given by Fc — ^ij^i^ji where 

(i, j) is a pair of any layers and V{xi — Xj) are Coulomb 
energy between two layers if z 7^ j and band energy if 




FIG. 2: The doping concentration near the surface. The 
average doping concentration inside the bulk is chosen to be 
0.2. 



i = j- Since the physics is very general, we simply assume 
the sheet-charge model,— in which the uniform distribu- 
tion of carriers in the Cu-O plane forms two-dimensional 
sheets of charge. The charge in the A-0 layers, except 
the one exposed on the surface, are fixed and determined 
by doping concentration. The charge in the layer at the 
surface can be freely adjusted since it is exposed to out- 
side. Therefore, we can consider a uniformly doped su- 
perconductor with negative charge accumulated on the 
surface with density —Sq. This surface charge induces a 
non-uniform doping distribution of charges Si in the ith 
Cu-O plane. Then the band energy under the carrier 
distribution, pi = Si/ 5 along the c axis can be written as 



(1) 



with m* w Anie And the electric field is E{k) — 
(47r/e)(5(l — Pj); where e is the background dielec- 
tric constant. The relevant electrostatic energy density 



is Ueii) = (e^ 



The total charge carrier en- 



ergy can be written as 



(2) 



where C = Am* e'^ do / fi^ e and do is the distance between 
the adjacent planes. Using the values do — 11.69A and 
e « 12eo, one obtains C « 7.4. 

In the presence of a competing order, we adopt the 
Ginzburg-Landau free energy of a multi-layer system 
given in Reii^, 



j 



[aiS,)\^,\^ + X\^,\^-p,i^,r,+i+c.c.) 
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(3) 



where a{Sj) and a'{Sj) are functions of the doping Sj for 
each layer and all other parameters are assumed to be 



3 




1 2 3 4 5 

i 



FIG. 3: The configurations of (a) tlie competing and (b) su- 
perconducting order parameters wlien tlie voltage is applied. 
V is the voltage dropped near the surface. The curves are cal- 
culated at V = 100,0, —100 meV respectively. The absolute 
value of V is expected to be less than the bias voltage of the 
sample with respect to the tip. 



constants. (f>i and ipi are the competing order and su- 
perconducting order parameters on the ith Cu-0 plane 
respectively. The total free energy of the system is given 
by Ft — Ff. + Fq. For a multi-layer system, in general, 
we have to resort to numerical method to minimize the 
total free energy. However, as an approximation, one 
can minimize Fc first with respect to Si and obtain the 
equilibrium charge distribution {6i}. After obtaining the 
charge distribution by minimizing Fc, we can minimize 
Fo with respect to 4'j a-nd 4>j to find the configurations of 
the competing and superconducting order parameter as 
a function of layers. This simple approach is proved to be 
reliable since Coulomb energy is usually larger and only 
a few layer closest to surface are effected which we will 
show later. In fact, we have performed an independent 
numerical calculations to verify this by using of a varia- 
tional Monte Carlo calculation to search the global min- 
imum of the total free energy. The quantitative results 
for the distribution of charge density and order parame- 
ter configuration are only slightly different from the ones 
obtained through the above approximated approach. 

Without external voltage, it is obvious that the A-0 
layer at the surface prefers zero charge density, i.e. 6o = 
0, in order to minimize the charge energy. Therefore, 
the first Cu-0 near the surface will be underdoped in 
general. In the presence of external voltage, additional 
electric field can be built near the surface, which leads to 
an effective modification of the charge at the surface, i.e., 
S{V) = 5o + f{V), where f{V) is zero iiV = and carries 
the same sign as the voltage V. Our natural assumption 



here is that 

f{V) = coV. (4) 

If we assume that the entire voltage is dropped in the 
competing order region which is roughly in one or two 
layers along c-axis, the value of cq can be estimated to 
be around ^ which is 0.09{V^^) for the lattice constant 

a = 4A, do = 11.69A, and e = l2eo- In the underdoped 
superconductor {5j < 0.2, for all j), wc choose a{Sj) — 
10{Sj - 0.3), a'{Sj) = 27{Sj - 0.22), A = A' = 1 and 
g = 1.2. This choice of the parameters gives the right 
shape and the magnitude of the superconducting dome 
for the generic superconducting phase diagram^ 

III. TUNNELING SPECTRUM 

The doping concentration at zero bias near the surface 
are calculated and shown in Figure [21 where the average 
doping concentration in the bulk is taken to be 0.2. This 
figure tells that only the first three layers closest to the 
surface are effected. The doping concentration recovers 
very quickly to the average values inside the bulk. The 
result also simplifies the numerical calculation to search 
a global minimum for the total free energy since only 
variables in a few layers need to be considered. 

Figure 121 shows how the configurations of the order pa- 
rameters respond to the applied voltage V, which is the 
voltage dropped near the surface. The curves in the fig- 
ure are calculated a.t V — 100, 0, —100 meV respectively. 
The absolute value of V is expected to be less than the 
bias voltage of the sample with respect to the tip. The 
positive V corresponds to electrons tunneling into the 
sample. As the applied voltage increases from negative 
to positive, the competing order on the second and third 
layers increases whereas superconducting order decreases 
on the second layer. In the deep inside the sample, the 
superconducting order is uniform. Near the surface, how- 
ever, the competing order parameter is dominant over the 
superconductivity. This means that the depletion region, 
where the competing order parameter (p is dominant, in- 
creases with the increasing voltage from negative to pos- 
itive. The charge depletion region can be defined by the 
average depth of the competing order parameter, l^, 

In Figure01 we plot the length of charge depletion region. 
Id, as the function of the voltage. The doping concentra- 
tion is chosen to be optimal. 

In tunneling experiments, the tunneling matrix ele- 
ments usually are assumed to be constant. However, the 
presence of the charge depletion region obviously leads to 
the variation of the matrix elements. The increase of the 
charge depletion region near the surface with the voltage 
bias V implies the decrease of the conductance. Since 
the matrix elements are determined by the overlapping 
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FIG. 4: The length of the charge depletion region as the func- 
tion of the voltage bias. The doping concentration at bulk is 
set to be optimal. 
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FIG. 5: The tunneling spectrum calculated by adding the 
effect of the depletion region in cuprates in superconductring 
states. The solid line in the inset shows the difference be- 
tween the spectrum in negative voltage and one in positive 
voltage. The dashed dots are extracted experimental values 
from Refi^. 



of wavefunctions between the tip and the measured ma- 
terial, the change of the competing order near surface can 
result in a significant variation of the tunneling matrix 
element. A slight variation of competing order near sur- 
face could lead to large change of the tunneling-matrix 
element. 

The effect is universal as soon as a competing order 
exists close to superconducting phase, regardless of the 
nature of the competing order. The effect is large if the 
competing phase is insulating. One promising compet- 
ing order is the antiferromagnetic (AFM) order, ^ which 
is an insulating phase. It has been suggested that the 
transition between AFM and SC orders are controlled 
by chemical potentialjiiSS which can be changed by the 
doping concentration. Moreover, it has also been ar- 
gued that the superconducting materials are locked at 
the AFM and SC transition point in the large range of 
doping concentration)2£ which inevitably results in a mi- 



croscopic phase separation. In this case, near the surface, 
the material should be very sensitive to external electric 
field since it effectively changes the chemical potential. If 
the competing order is metallic, such as d-density wave 
(DDW),^ the effect can be small. The value of asymme- 
try is also expected to depend on the disorder which is 
generically present in the cuprates. Although it is hard 
to make an quantitative estimation since the different 
competing orders will lead to different magnitude of the 
asymmetry, we can make a reasonable assumption to ob- 
tain quantitative result to see whether the experimental 
result can be naturally explained. 

In a simple approximation, we can assume that the de- 
pletion region plays as a tunnelling barrier in the STM 
measurement between the STM tip and the supercon- 
ducting phase. The energy of the barrier is determined 
by the single-particle energy gap Eg . Since the depletion 
region is close to half filling, Eg should be expected in 
order of U in the microscopic Hubbard model. This ap- 
proximation leads to an exponential dependence of the 
tunneling matrix element as a function of the length of 
the depletion region. 



where Go{V) is the tunneling conductance for a d-BCS 
state which is assumed to be symmetric for negative and 
positive bias voltages i.e. Go{V) = Go{—V), do is the 
lattice constant along c-axis and ld{V) is the depletion 
length in the unit of lattice constant in c axis. In Figure^ 
we calculated a typical spectrum following the above ap- 
proximation. The asymmetry between positive and neg- 
ative bias voltage, G{—V) — G{V) is also shown in the 
inset. Amazingly, the experimental data can be fitted to 
this crude approximation. The dotted points in the in- 
set of figure[Slare the extracted experimental values from 
Refii4 Using the result of ld{V) calculated by Equation 
101, there is only one fitting parameter Eg. For the ex- 
perimental resultji^ we obtain Eg = 0.79 eV, which is 
indeed in the order of U. Therefore, this result shows 
that the suppressed superconducting order parameter by 
the dominance of competing order parameter close to the 
surface of the sample can explain the asymmetry of the 
tunneling conductance background of high temperature 
superconductors. 



IV. CONCLUSION AND DISCUSSION 

This theory provides several explicit predictions. First, 
there is an onset temperature for the tunneling asymme- 
try since there is onset temperature for the competing 
order. In cuprates, a natural onset temperature could be 
T* associated with the pseudogap. We expect that the 
tunneling asymmetry disappears above T* . Secondly, the 
tunneling asymmetry should be larger if the bias volt- 
age is larger. It is a monotonically increasing function 
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of the voltage, which is consistent with the experimen- 
tal result of the tunneling asymmetry as the function of 
voltage. The experimental result shows that the asym- 
metry is small in low voltage and becomes larger and 
larger as the voltage becomes largeJ^ The theory predicts 
that there could be a large enhancement of the asymme- 
try at high bias voltage if the nature of the transition 
between the competing order and superconductivity is 
the first order transition. For example, the transition 
between AFM and SC is typically first order. Finally 
the asymmetry which has recently been observed in the 
point-contact spectroscopy of the heavy fermion super- 
conductor CeCoIngi^ can also be naturally explained by 
this effect since the material in this experiment is close to 
the boundary of the AFM and SC phase transition. Our 
theory predicts that the asymmetry should be weakened 
when the material moves away from the boundary toward 



superconducting zone by tuning the pressure. 

In conclusion, we have presented that in the presence 
of competing orders at surface, the tunneling matrix el- 
ements can strongly depend on externally applied volt- 
age. This effect can lead to a natural explanation of the 
tunneling asymmetry observed in the STM experiments. 
We have shown that this effect could exist universally in 
materials where the superconducting phase is close to a 
competing order and our prediction can be easily tested 
experimentally. 
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